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Abstract: The giant radio galaxy M87 is the only extragalactic non-blazar object which has been detected 
as a source of very high energy 7-rays. It represents a unique opportunity to study the phenomena of 7- 
ray emission from a nearby AGN. In this paper we report preliminary results from the observations of 
M87 taken with the imaging atmospheric Cherenkov telescope array VERITAS in February, March and 
April 2007. An excess of photons above an energy threshold of 250 GeV is measured with a statistical 
significance of more than five standard deviations. 



Introduction 

M87 is a nearby giant elliptical galaxy (~16 Mpc) 
which lies near the center of the Virgo galaxy 
cluster. It is a powerful radio source (Virgo A), 
classified as a Fanaroff-Riley class I (FR I) radio 
galaxy [9]. Its core is an AGN powered by a 
supermassive black hole of about 3.2 x 10 9 M Q 
[14], emitting the first-detected plasma jet [7] 
which extends over several kiloparsecs. According 
to the unified scheme [18] of BL Lacertae objects 
(BL Lac) and FR I radio galaxies, M87 should be 
of the same astrophysical nature as BL Lac but 
with its jet not pointing along the line of sight. 
The critical sight angle between BL Lac and 
FR I galaxies is not well known. Some models 
propose an M87 jet misalignement around 30° [3] 
but superluminal motion observed by HST [5] in 
the M87 jet suggest an orientation within 19° of 
the line of sight. The TeV emission from M87 
establishes an important link between M87 and 
BL Lac objects, which may have implications for 
the unified scheme, since all other extragalactic 
sources detected in the TeV regime are BL Lac 
objects [11]. 

M87 has been observed over a broad range of ener- 
gies from radio waves to 7-rays. Its jet is resolved 
in radio, optical and X-ray regimes and shows sim- 
ilar morphologies at all wavelengths. This emis- 



sion is understood as the synchrotron radiation 
from high energy electrons in the jet. The same 
electron population could also be responsible for 
the TeV 7-ray emission through the inverse Comp- 
ton process [17]. Models involving protons [16] 
can also explain the TeV emission. 

In fact, M87 is also considered as a possible 
source of ultra-high-energy (~10 20 eV) cosmic 
rays [4]. As tracers of high-energy hadrons, 
very-high-energy 7-rays could provide important 
clues to help determine the origin of extragalactic 
cosmic rays. 

The first detection of M87 in the TeV regime was 
reported by the HEGRA collaboration [1] with a 
statistical significance of 4.1 standard deviations. 
This was confirmed by the HESS collaboration [2] 
which also reported rapid variability and an unex- 
pectedly hard spectrum. The imaging atmospheric 
Cherenkov technique provides insufficient angular 
resolution (~ 3 arcmin) to resolve the M87 emit- 
ting region (core, jet...) but the day scale variabil- 
ity measured by HESS suggesting a very small re- 
gion, most likely close to the core. The observa- 
tion of superliminal motion of extremely compact 
radio structures within the HST-1 region has how- 
ever recently provided indications the TeV emis- 
sion could originate from that region, 0.82" (60 pc 
projected) from the core [6]. 
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Observation with VERITAS 

VERITAS is an array of four imaging atmospheric 
Cherenkov telescopes with 12 m-diameter mirror 
and 3.5° field of view camera. Telescopes details 
can be found in these proceedings [15]. The array 
is located in southern Arizona where M87 reaches 
more than 70° in elevation. 

The observations in the direction of M87 were car- 
ried out during 51 hours from February to April 
2007 at elevations from 55° to 71°. All data were 
taken in "Wobble" mode, tracking M87 with a 0.5° 
offset (successively north, south, east, west). VER- 
ITAS was still under construction and most of the 
observations were performed with a partial array of 
three telescopes. In order to check the source vari- 
ability, we processed all the data in the same way, 
using only the first three telescopes even when the 
fourth telescope was in use (6% of the data set). 
We selected usable data on the basis of the weather 
conditions and the raw trigger rate stability. About 
90% of the data pass our quality selection cuts (cf. 
Table 1). 

Data Analysis 

The data have been analyzed using independent 
analysis packages (see [8] for details on the anal- 
ysis). All of these analysis yield consistent re- 
sults. The presented results are obtained with the 
GrISU package [10]. This analysis chain differs 
from those described in [8] by fixed thresholds in 
the image cleaning (same thresholds for all pixels) 
and by the background-rejection cuts. In this anal- 
ysis, the background-rejection cuts were optimized 
for the detection of a point-like source with a lu- 
minosity of 10% of the Crab Nebula. The signal 
from the source ("ON") is the number of events 
with a reconstructed direction less than 0. 14° from 
the source direction. The backround rate ("OFF") 
is estimated using seven "mirror" regions which 
have the same area and the same distance from the 
camera center (0.5°) as the "ON" region (reflected 
region method [8]). The statistical significance is 
calculated using the Li and Ma formula [13] (with 
a = 1/7). 

This analysis, applied to the Crab Nebula data 
taken with the 3-telescope-array configuration at 
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Figure 1: Excess significance map of the M87 sur- 
rounding region. 



~78° elevation in February 2007, results in an 
excess of 5.9 7/min for a background rate of 
0.75 evt/min (Li & Ma excess significance of 
27.5 a/Vh). In this preliminary paper, we express 
the measured 7-ray flux in Crab units. Simulations 
and data show the excess rate does not vary signif- 
icantly (< 10%) between 78° and the M87 obser- 
vation elevation. 

Preliminary results 

VERITAS detected an excess of 263 events from 
the direction of M87, corresponding to a statistical 
significance of 5.1 standard deviations (secondary 
analysis yield also a significance above 5 standard 
deviations). The energy threshold of our analy- 
sis at the M87 observation elevation has been es- 
timated with Monte Carlo simulations at 250 GeV 
The time averaged excess rate corresponds to a 7- 
ray flux of 1.7% of the Crab nebula flux. 

Figure 1 shows the map of the excess signifi- 
cance in the surrounding region of M87. The ex- 
cess position is: R.A. 12 h 30™ 47.5 s ± 11 s , Dec. 
+12°23'23" ± 1'48", consistent with the position 
of M87 (R.A. 12' 1 30™ 49.4 s , Dec. +12° 23' 28"). 
Figure 2 shows the distribution of the square of the 
angle 8 between M87 and the reconstructed direc- 
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Date 


Observation time 


Results 


raw 


selected 


excess rate 


significance 


12-25 Febr. 


17.0 h 


15.5 h 


0.14± 0.04/min 


4.2cr 


12-20 March 


23.5 h 


20.4 h 


0.07 ± 0.03 /min 


2.6 a 


07-19 April 


10.7 h 


8.3 h 


0.10 ± 0.05 /min 


2.2 d 


TOTAL 


51.2h 


44.2 h 


0. 10 ± 0.02 /min 


5.1 a 



Table 1: Data set toward M87 taken with VERITAS in spring 2007 
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Figure 2: Crosses show the 9 2 distribution from 
M87 and the solid line from off-source directions. 
The dashed line is the 8 2 distribution expected for 
a point-like source with a backround level repre- 
sented by the dotted line. 
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Figure 3: Nightly average (crosses) and global av- 
erage (dashed line) of the M87 7-ray flux above 
250 GeV as a function of the observation day of 
year. 
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Figure 4: M87 flux reported by very high energy 
7-ray telescopes as a function of observation year. 
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tion. The shape of the excess is compatible with a 
point-like source (dashed line). The upper limit to 
the source extension has a 2.3 arcmin radius. 

Figure 3 shows the M87 light curve night-by-night 
during the three months of observation. No signif- 
icant variability is observed. The constant flux fit 
X 2 per free parameters is: 21.4/22. The maximal 
deviation of the rate from the average (February 21 
and 22) reaches only 3.4 standard deviations. We 
cannot yet confirm the variability on a 2-day scale 
reported by the HESS collaboration. 

Figure 4 shows the M87 flux recorded during the 
last 10 years by HEGRA [1], Whipple 10 m [12], 
HESS [2] and VERITAS expressed in Crab units. 
The 2007 flux reported in this paper is at an average 
level. 

Summary and Conclusion 

VERITAS confirms the 7-ray emission from M87 
above 250 GeV The measured flux is below 2% 
of the Crab Nebula flux, comparable to earlier ob- 
servations with HESS and HEGRA. This result 
demonstrates the capability of VERITAS for de- 
tecting a faint source in 50h of observations. 
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